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The Fast-Sup erf ast Transition in the Sleeping Lagrange Top 

Tudor S. Ratiu* Miguel Rodriguez-Olmos | Miguel Teixido-Roman * 


Abstract 

The fast-superfast transition is a particular movement of eigenvalues found in [4] when studying 
the family of sleeping equilibria in the Lagrange top. Although this behaviour of eigenvalues typically 
suggests a change in stability or a bifurcation, in this case there is no particular change in the 
qualitative dynamical properties of the system. Using modern methods, based on the singularities 
of symmetry groups for Hamiltonian systems, we clarify the appearance of this transition. 


MSC 2010: 70H05; 70H14; 37J20 

1 Introduction 

The Lagrange top is one of the most well-known simple mechanical systems with symmetries. It 
consists of an axisymmetric rigid body with a fixed point moving under the influence of a constant 
gravitational force. One of the important characteristics of this system is that it has two symmetry 
groups: the rotations around the axis of gravity and the spin about the axis of symmetry of the 
top. In fact, the energy and the conserved quantities associated with these symmetries make the 
Lagrange top a Liouville integrable system (see, e.g., ®)- 

A special configuration of the system is the sleeping top, when the top is pointing upwards, so that 
the axis of gravity and the symmetry axis coincide. This position has non-trivial isotropy, since there 
is a continuous group of symmetries that leaves this configuration invariant. This fact is ultimately 
responsible for many of the dynamical properties of a sleeping top. In this work, we focus on this 
particular configuration, but with all possible angular velocities. These solutions are examples of 
relative equilibria, motions of the system which evolve within the symmetry direction. Generally 
speaking, relative equilibria act as organising centres of the dynamics of a Hamiltonian system with 
symmetries and, therefore, their study gives important information about the qualitative behaviour 
of the dynamical flow which, typically, is impossible to obtain analytically. 

Unlike for generic flows, relative equilibria of Hamiltonian systems can be characterised as critical 
points of a certain function depending on a parameter. This approach is very useful, because it 
gives necessary conditions for bifurcations of branches of relative equilibria, based on singularity 
and critical point theory. Indeed, when the second variation of this parameter dependent function 
becomes degenerate, a new branch of relative equilibria may bifurcate. Typically, one expects that 
at the bifurcation point, a transfer of stability will occur from the original to the bifurcating branch. 
However, in the Lagrange top, due to the existence of continuous isotropy, this is not at all the 
case. For a large part of the family of sleeping Lagrange tops, every point is a bifurcation point to 
precessing solutions, without the original branch loosing stability. The possibility of a bifurcation 
from a stable branch to a branch that is also stable was already pointed out by Cartan [1J. 

These and other phenomena are well known and have been studied in depth from a geometric 
perspective in [2| , where the stability range of sleeping Lagrange tops and their possible bifurcations 
are computed and the spectral analysis of the linearised dynamics along the sleeping Lagrange top 
family is carried out. Assuming an oblate body, i.e., the largest moment of inertia is along the 
symmetry axis, in [3], a diagram similar to the one showed in Figure [l] is presented. 
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Figure 1: Eigenvalues of the sleeping top as a function of the angular velocity A in an 


oblate top 


When the angular velocity A is zero, the eigenvalues are a real pair. As A is increased, the 
eigenvalues form a quadruple and when the angular velocity attains the fast-slow critical value 
t/_ s , the eigenvalues move onto the imaginary axis and the system becomes stable. This fast-slow 
transition corresponds to a Hamiltonian-Hopf bifurcation (see [3]). The surprising fact noticed in 
3] is that if the angular velocity is further increased, two of the eigenvalues cross at zero moving 
along the imaginary axis for A 2 = r/_ s /. This point was called the fast-superfast transition. Usually, 
such an eigenvalue crossing implies the existence of a bifurcation or a change in stability, but from 
a dynamical point of view, this fast-superfast point has no special behaviour with respect to nearby 
points in the sleeping top family. On the other hand, one could argue that since all the points in 
the stable regime exhibit a bifurcation to the precessing branch, a zero eigenvalue at each point in 
the stable range is expected. However, only a zero eigenvalue is observed in the linearisation at the 
fast-superfast transition point. 

If the body is prolate, i.e., the shortest moment of inertia is along the symmetry axis, the evolution 
of the eigenvalues is shown in Figure[2] The motions of eigenvalues shown in both figures correspond 
to those obtained in Figure 3 of jj. 
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Figure 2: Eigenvalues of the sleeping top as a function of the angular velocity A in prolate top 


As in the oblate case, the fast-slow transition is a Hamiltonian-Hopf bifurcation between the 
stable and unstable regimes, but note that in this case, no fast-superfast transition occurs. It can be 
shown that the qualitative dynamical behaviour of the system does not depend on the oblateness or 
prolateness of the body, so the linearisation should not behave differently in the oblate and prolate 
cases. In particular, the existence of a fast-superfast transition at a unique point in the oblate but not 
the prolate case, questions our understanding of the qualitative dynamics of the sleeping Lagrange 
top family. 

In this paper, we show that this fast-superfast transition is an artefact due to the existence of 
continuous isotropy in the sleeping Lagrange top. This has as a consequence that the linearisation 
done in [3 has implicitly chosen one angular velocity representative among an infinite number of 
possibilities for the equilibria under study. In fact, we show that for each stable sleeping equilibrium, 
we can choose the angular velocity in such a way that the linearisation has a double zero eigenvalue 
crossing and therefore exhibits a fast-superfast transition, for both the oblate and prolate cases. 


2 The Setup 

The Lagrange top is the symmetric Hamiltonian system defined by the tuple 

(T*SO(3), w c , T 2 , J, h), 

where the phase space is the cotangent bundle T*SO(3) of the proper rotation group SO(3), equipped 
with its canonical symplectic form ui c - If we use the identification, given by the right translation 

V := SO(3) xl 3 -x T*SO(3) 

(A, 7T ) —X (A, 7rA), J- el 3 , (1) 
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then the toral symmetry group T 2 acts on on V as 

((0i,02),(A,tt)) ^ (e eiS5 Ae- e2g3 ,e fllg5 7r), (6>i ,8 2 ) G S 1 x S 1 . 

The associated canonical momentum map for this action is given by 

J(A, 7r) = (n ■ e 3 , — 7r • Aes) G R 2 . (3) 

Finally, the Hamiltonian function h : V —► R is 

h( A, 7r) := mgle 3 • Ae3 + ^7r ■ Ia 1 ?]-, 

where Ia := AIA -1 and I := diag(/i, Ji, I 3 ) is the inertia tensor of the rigid body in the reference 
configuration with respect to a principal axes body frame. Notice that, in this reference configuration, 
the body symmetry axis is e 3 . Throughout this article we are using both R 3 and the space so (3) 
of 3 x 3 antisymmetric matrices to represent the Lie algebra of SO (3), identified via the usual Lie 
algebra isomorphism given by the hat map 


“ : R 3 

—>• so(3) 

0 

x 3 

X2 

X = (X1,X2,X 3 ) 

X = 

X 3 

0 

—XX 



-X2 

Xx 
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In right trivialisation we have <5A = 59A G 7aSO(3). Therefore, a tangent vector v G T(a >7[ \V is 
represented as (58, Sir) G R 3 xR 3 . 

From this, it is easy to check that the expression for the symplectic form ujonV (the pull back 
of <jj c to V given by 0 ) is given by 

co(A, 7r)( (58i, 5m) (59i, 5ty2)) = 5 n 2 ■ 50i — Styx ■ 582 — 7r ■ ( 58i x 682 ). (4) 


3 The sleeping Lagrange top 

In this model, the sleeping Lagrange top is a relative equilibrium at the phase space point z = 

(A, 7 r) := (Id, XI 3 e 3 ) with arbitrary velocity (£l,£r) G R 2 = LieTF 2 satisfying — ^r = A. This 
makes the set of all admissible velocities at each relative equilibrium a one-parameter family. In 
order to see this, and using standard arguments from geometric mechanics (see, e.g., [5], | 6 ]), we 
have to check that the augmented Hamiltonian 

7r ) : = mgle 3 ■ Ae 3 + - n ■ e 3 ^ L + n ■ Ae 3 £ R (5) 

has a critical point at (Id, XI 3 e 3 ). The derivative of h^ L ^ R ) at an arbitrary point (A, n) £ V is 

d/i( eti 5 R )(A, n)(59, 5n) = mgle 3 -(59x Ae 3 )+j7r-lA 1 7r-|-7r-((56lxlA 1 7r)— ^L5n-e 3 +^R5n-Ae 3 +^ R n-(59x Ae 3 ) 

At (A, 7 r) = (Id, XI 3 e 3 ) this becomes 

dh (€L,£ r )( Id > A7 3 e 3 ) (58, Sty) = (A - (£ L - £r))5ty ■ e 3 ( 6 ) 

which vanishes precisely when — (,r = X. Therefore, according to (|2|, the dynamical evolution of 
z = (Id, A I 3 e 3 ) is given by z(t) = (e tXe3 , XI 3 e 3 ). 

In order to study the linearisation of the Hamiltonian system at the sleeping equilibrium, we need 
to compute several more geometric objects. First, from |3|, it is clear that 

M:= J (z) = (A/s,-A/s). 

Second, since T 2 is Abelian, the coadjoint stabiliser of g is = T 2 . Third, from Q, we find that 
the stabiliser of the phase space point 2 = (Id, XI 3 e 3 ) is 

G z = S lD := {(8, 8) G T 2 : 8 G S' 1 }. 

A normalised basis for its Lie algebra is Aj(l, 1), i.e., 

= span 

Define m := span |^j(l, —1)| which is a complement to 0 Z in 0 , i.e., 

0 = m © 0 .. (7) 
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According to this direct sum decompostion, the velocity of the sleeping Lagrange top takes the form 

(&.&) = + »?(!. 1 ). ( 8 ) 

where r/ £ R is arbitrary. Since G z is a continuous subgroup of positive dimension, 0 reflects 
the fact that the velocity of the relative equilibrium z, for which the first variation § vanishes, is 
defined only up to an element of g 2 . However, its projection onto the subspace m, according to 
the splitting 0, sometimes called the orthogonal velocity of the equilibrium, is unique. 


4 The Symplectic Slice 

The linearisation of the Hamiltonian system at the relative equilibrium a = (Id, A/ 3 e 3 ) is given by 

L := n^d 2 h^ L , (R) (ld,XI 3 e 3 )\ N € sp(N), (9) 

where ( N, Qn) is a symplectic vector subspace of T Z V defined by an arbitrary G z -invariant splitting 

kerT^J = T Z (G^ ■ z) © N, 

fljv = uj(z)\ n , and d 2 h^ L ^ R ' ) (Id, AJ 3 e 3 ) denotes the Hessian of h^ Lt ^ R ) at the equilibrium point 
(Id, A/ 3 e 3 ). The vector space (N, fljv) is often called a symplectic slice. 

In order to compute the symplectic slice, we start by studying the fundamental vector fields 
(£i>£ 2 )p for the action 0- which are given, in the right trivialisation 0,by 

(£i,£2)t’(A,7t) = (£ie 3 — £ 2 Ae 3 ,£ie 3 x ir). 

So, at the sleeping equilibrium (Id, AJ 3 e 3 ), we have 


T(id,AJ 3 e 3 ) (G> • z) = span {(e 3 , 0)} , 

where p = J(z) = (A I3, —A I3) and the coajoint isotropy subgroup is G M = T 2 (see 0. 
Since the derivative of the momentum map is 

T( a,tt) J(<50, Sit) = (5n ■ e 3 , — 5n ■ Ae 3 — 7r • ( Sd x Ae 3 )), 

at the sleeping Lagrange top point (Id,A/ 3 e 3 ) this becomes 

T(id,A/ 3 e 3 ) 3{50,5n) = (Stt ■ e 3 , -Sn ■ e 3 - A J 3 e 3 ■ (60 x e 3 )) 


and, therefore, 


kerT (Id ,Ar 3 e 3 )J = {{80, Sn) : Sn _L e 3 }. 
Hence, a possible choice for the symplectic slice at z is 


N = {(60, Sir) : 60, Sir ± e 3 }, 


with symplectic form 

Qn ((60 1 , (Itti), (802,8712)) = 87V2 ■ 80 1 — Sir 1 • 802 — AJ 3 e 3 • (50 1 x 80 2 ). (10) 

It can easily be checked that N is indeed G z -invariant and that the isotropy group G z = S lD acts 
on N by 

ip ■ (80, Sir) = (e iv,S5 60, e^Sir), p £ S 1 ^ S lD . 


5 Stability 

In this section, we briefly reproduce the classic stability result for a sleeping equilibrium of the 
Lagrange top in a formulation adapted to the general framework of this article. This is necessary for 
the subsequent bifurcation analysis relevant to the fast-superfast transition. 

The second variation of the augmented Hamiltonian 0 is 

d 2 h(£ L , £r)((80i, fori), (802,Stt 2 )) = 

+ 

+ 


(-mgl - \ 2 13 + A 2 -2- - A h(, R )(80i x e 3 ) ■ (S0 2 x e 3 ) 

Ji 


A 1- 


+ £,r [57Ti • (80 2 x e 3 ) + Sir 2 ■ (80 1 x e 3 )] 


Sir 1 ■ I 1 5 ir 2 - 
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Let 


( 11 ) 


ui = (ei,0), u 2 = (e 2 ,0), u :i = (0,ei), u 4 = (0,e 2 ) 
be a basis of N. Using it and (|8j, we get = g — A/2. Thus, at the sleeping top equilibrium, 


d 'fys_L,£n)( Id > A/ 3 e 3)| Ar = 


A 0 OB 
0 A -B 0 
0 —B C 0 
BO 0 C 


( 12 ) 


where 


A:m-mgl+^(2I 3 -Ii)-\I 3 r,, B := ±(2I 3 - Ii) - g, C := j-. 


The eigenvalues of this matrix are 

1 


<T± := \ ((A + C) ± y/(A + C) 2 - 4 [AC - B 2 )) . 


Thus, o± > 0 <=> AC - B 2 > 0. 

Recall from the theory of stability of Hamiltonian relative equilibria (see [7] for a reference 
appropriate to the isotropy-based approach consistent with this article), that in order to guarantee 
nonlinear stability it is enough to find an admissible velocity such that ( |12| ) is definite. For a sleeping 
top equilibrium of the form z = (Id, XI 3 e 3 ), this is equivalent to saying that it is stable if we can find 
g £ R such that AC — B 2 > 0. 

Note that 

AC-B a = V + - 

il 2ii 1 i 4 

Therefore, for a Hxed value of A, this is a quadratic function of r/ that attains its maximum value 
when rj is equal to 

» h — h . 

71 '= ~ 2 ir X - 

Substituting this value of r/ in the previous expression yields 

A 2 /f — Amglh 


AC -B 2 


47J 


and we conclude that if 


A > 


Amglli 


T 2 

1 3 


then the sleeping top equilibrium (Id, XI 3 e 3 ) £ V is nonlinearly stable. This is the classical fast top 
condition. Notice that if 

a 2 Amglli 


l 2 

then the stability test is inconclusive and this regime must be studied by linearisation methods, as 
we do in the next section. 


6 Linearisation 

Formula (|10| easily implies that the expressions of the linear symplectic form S7 jv and its inverse, 


relative to the basis 1111, are 



' 0 

-XI 3 

1 

o' 


'0 

0 

-1 

0 

fljv = 

A I 3 
-1 

0 

0 

0 

0 

1 

0 

Q —1 — 

» ^ l N — 

0 

1 

0 

0 

0 

0 

-1 
—A I 3 


0 

-1 


0 


0 

1 

A I 3 

0 


so the matrix of the linearised system is 


L = 


0 

—B 

A - XI 3 B 


B 

0 

0 


-c 

0 

0 


0 

-C 

B - XI 3 C 

0 


0 A — XI 3 B -B + XI 3 C 
The characteristic polynomial of L is 

t 4 + (2 {AC - B 2 ) + (2 B - CI 3 X) 2 ) t 2 + (AC - B 2 ) 2 , 


(13) 
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which, after some manipulations, can be written as 


S + E ' 


2 II 

where 

E :=/ 3 A-/i( 27 ? + A), 
Thus, the eigenvalues of L are 


+ F t 2 + ( E2 ~ F ) 2 


16/4 


F := 7|A 2 - Amglh. 


e±,± ■■= ^(±e±Vf). 


( 14 ) 


Notice that for any A, that is, for any sleeping top equilibrium, an admissible velocity (equiva¬ 
lently, a value of ?;) can be chosen so that E takes any real value. Also, F is positive if and only if 


A corresponds to an equilibrium for which there is an g making the matrix (12 I definite. Therefore, 
we conclude that: 

• For any A 2 < 4m 7 iJl , the eigenvalues of the linearisation have non-zero real part and therefore 

J 3 

the sleeping equilibrium is unstable. E can be chosen so that the imaginary part is zero (real 
double pair) or non-zero (complex quadruple) at each point of the unstable regime. 

• For any A 2 > 4m A ;/l , the eigenvalues of the linearisation are purely imaginary. At each point of 

J 3 

the stable range, E can be chosen so that we have 4 distinct non-zero eigenvalues, two double 
imaginary non-zero eigenvalues, or an imaginary pair and a double zero. 

Notice that the last property guarantees the existence of one imaginary pair plus a doble zero 
of the linearised system at each point of the stable range, provided a suitable admissible velocity at 
that point is chosen. This shows that the fast-superfast transition actually happens at each point of 
the stable regime if we use the freedom in the isotropy Lie algebra when linearising the system at a 
sleeping equilibrium. 

We now study the relationship between our approach and the results in [4]. First, we notice 
that in order to obtain a zero crossing of eigenvalues along the sleeping top equilibrium family (the 
fast-superfast transition) the condition to be satisfied is 


which is equivalent to 


E 2 = F, 


(2/ 3 — Ji)A 2 + 4(7 3 — Ii)r/X — 4/i g 2 — Amgl = 0. 


(15) 


Notice that the moments of inertia inequalities 7, < Ij + A for any distinct imply that this 

quadric is a hyperbola in the (A, ??)-plane for any choice of inertia tensor. As we have seen before, 
one can choose g such that (A, g) lies on the hyperbola given by (151 if and only if F > 0 (stable 
sleeping equilibrium), which holds for both the prolate and oblate cases, as opposed to was observed 
in g]. 

A reason for this apparent disagreement is the following. The approach taken in 2], and based on 
reduction of the right S' 1 -action, corresponds, in our setup, to making the permanent choice g = A/2 
along the family of sleeping equilibria. Therefore, in 2], along this family, 

(£l,£r) = (A, 0). 


Using this choice of velocities, we conclude that (151 is equivalent to 


mnl 

h - 7, ’ 


(16) 


which corresponds, in the axisymmetric case Ii = I 2 , to the fast-superfast conditions (4.24) of 
2]- Notice that, in the (A, 77 )-plane, for a prolate top I 3 <7 1 , the line g = A/2 never meets the 
hyperbola (151. However, for an oblate top J 3 > h, the fast-superfast transition is found at the 


unique point (161, which is the intersection of this line with the hyperbola. Both observations are in 


total agreement with 2]. 

These considerations are illustrated in the following figures. In Figure[3] the different possibilities 
for the eigenvalues of the linearised system are shown in the oblate case. The line L \, corresponding 
to the choice g = A/2 of (4/ intersects the hyperbola E 2 = F precisely at the value A = r/_ s /. Notice, 
however, that it is possible to find a different relationship between A and g such that this intersection 
happens at any value of A greater than t/_ s , that is, the fast-superfast transition happens at each 
point of the stable regime in the family of sleeping equilibria for oblate tops. 

In Figure[4]the analogous situation for prolate tops is presented. Notice that in this case, the line 
Li corresponding to the velocity choice in 2] does not intersect the hyperbola E 2 = F. However, as 
is apparenr from the figure, and exactly as for the oblate case, a different linear relationship between 
A and g can be chosen so that this intersection exists for any value of A in the stable regime. 

In both figures we have shown the linear relationship between A and g corresponding to the 
condition E = 0. Along that line, the eigenvalues of the linearisation always come in pairs; the 
fast-superfast and the fast-slow transitions occur at the same point. 
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Figure 3: Oblate case 



Figure 4: Prolate case 


7 Conclusion 


The existence of continuous isotropy, and therefore an isotropy Lie algebra of positive dimension, 
at a relative equilibrium reflects the fact that the linearisation of the Hamiltonian vector field at the 
equilibrium, given by is not uniquely defined. The effect of this is that the set of all linearised 
vector fields should be considered in order to avoid missing pieces of information about the qualitative 
properties of the Hamiltonian flow near the equilibrium point. In the case of the sleeping Lagrange 
top, we have shown that the systematic study of the j),-parametrised family of linearisations given 


by (131 guarantees that the fast-superfast transition occurs for both prolate and oblate Lagrange 


tops. Furthermore, in both cases, this transition can be observed for all points in the stable range 
along the family of sleeping equilibria. 

It can be shown that the fast-superfast transitions are closely related to the bifurcations from 
sleeping to precessing relative equilibria, which are known to happen precisely at each point of the 
stable range of sleeping tops. This will be studied in detail in [S]. 
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